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Abstract 

Many cells contain non-centrosomal arrays of microtubules (MX), but the assembly, organisation 
and function of these arrays are poorly understood. We present the first theoretical model for the 
non-centrosomal MX cytoskeleton in Drosophila oocytes, in which bicoid and oskar mRNAs become 
localised to establish the anterior-posterior body axis. Constrained by experimental measurements, 
the model shows that a simple gradient of cortical MX nucleation is sufficient to reproduce the 
observed MX distribution, cytoplasmic flow patterns and localisation of oskar and naive bicoid 
mRNAs. Our simulations exclude a major role for cytoplasmic flows in localisation and reveal an 
organisation of the MX cytoskeleton that is more ordered than previously thought. Furthermore, 
modulating cortical MX nucleation induces a bifurcation in cytoskeletal organisation that accounts 
for the phenotypes of polarity mutants. Xhus, our three-dimensional model explains many features 
of the MX network and highlights the importance of differential cortical MX nucleation for axis 
formation. 
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INTRODUCTION 


Microtubules (MTs) are polar cytoskeletal filaments that can adopt different global net¬ 
work organisations to fulfil different functions. The vast majority of studies have focussed on 
understanding the architecture and function of MT arrays organised by centrosomes, such as 
radial arrays or the mitotic spindle. In contrast, much less is known about the organisation, 
assembly and function of non-centrosomal MT arrays despite their ubiquity in differentiated 
cell types, such as neurons, epithelia, hssion yeast and plants niEi. In a variety of cell-types, 
non-centrosomal MT arrays play an essential role in directing the subcellular localisation of 
mRNAs to spatio-temporally control gene expression. In neuronal dendrites, for example, 
MTs form bidirectional arrays with MT plus-ends both pointing away and towards the cell 
body Eli]. MTs and associated motor proteins were implicated in the transport of mR¬ 
NAs along dendrites Eli where their activity-dependent translation contributes to long term 
changes in synaptic function, neuronal circuitry and memory [6HH|. Similar overlapping, 
bidirectional MT arrays in the Xenopus oocyte have been proposed to mediate the trans¬ 
port of Vgl mRNA to the vegetal cortex, where it orchestrates germ layer patterning [a HO]. 
The Drosophila oocyte is probably the best-studied example of mRNA transport along non- 
centrosomal MTs. In this system, a diffuse gradient of MTs of mixed polarity is required 
for the localisation of bicoid and oskar mRNAs to opposite ends of the cell m- Despite the 
large amount of work, however, the organisation of the non-centrosomal MT cytoskeleton 
underlying this mRNA localisation is controversial [IMl], and its assembly and function 
are not understood. In stage 9 oocytes, MTs grow from most parts of the cell cortex into 
the volume da US] thereby giving rise to a complex, three-dimensional MT network without 
pronounced polarity along the anterior-posterior axis [lain]. In contrast to this apparent 
disordered organisation, bicoid and oskar mRNAs become reliably localised by Dynein and 
Kinesin to the anterior corners and to the posterior pole of the oocyte, respectively, thereby 
dehning the anterior-posterior axis [HdHHa]. 

The most pronounced feature of the MT cytoskeleton at stage 9 is a gradient of cortical 
MTs from the anterior to the posterior pole, where MT nucleation is suppressed by the 
polarity protein PAR-1 [^123] . Live imaging of oskar mRNAs and direct measurements of 
growing MTs showed that the MT network is mostly disordered with only a weak statistical 
bias of about 8% more plus ends pointing towards the posterior pole [H] . This bias vanishes 
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in absence of PAR-1 when MTs nucleate from all over the cortex [16]. While these hndings 
established a directionality of the MT meshwork for the hrst time, several questions remain 
unanswered. For example, mRNA transport and localisation is highly reproducible, raising 
doubts about whether the underlying cytoskeletal organisation can be mostly disordered. 
Moreover, Dynein-dependent transport to MT minus ends localises injected, so-called naive 
bicoid mRNA to the cortex closest to the injection site nzi, which is difficult to reconcile 
with a cytoskeleton that is simply biased towards the posterior everywhere. Similarly, the 
different behaviour of so-called conditioned bicoid mRNA, which localises specihcally to 
the anterior cortex irrespective of the injection site na, has remained unexplained. Finally, 
motor-driven cytoskeletal transport is not the only transport mechanism in oocytes. Kinesin 
moves only 13% of oskar mRNA at any given time, and the remaining 87% is subject to 
diffusion and slow cytoplasmic flows that are driven indirectly by Kinesin activity on the 
MT network [T3|. This raises the question if cytoskeletal transport alone is sufficient to 
account for mRNA localisation [2H [2S] and highlights the importance of distinguishing its 
contribution to localisation from the contribution of flows and diffusion [2B] . 

Here, we present the hrst theoretical model for stage 9 Drosophila oocytes. Based on 
the distribution of microtubule nucleation sites around the cortex, this model accurately 
reproduces the observed distribution of MTs and cytoplasmic hows in the oocyte. It reveals 
that the MT cytoskeleton is compartmentalised and more ordered than previously thought. 
By modelling the movement of mRNAs on this network, we also show that this MT organi¬ 
sation is sufficient to explain the localisation of oskar and naive bicoid mRNA. Finally, we 
show that modulation of MT nucleation gradients causes a bifurcation in cytoskeletal or¬ 
ganisation that explains mutant phenotypes. Thus, our results explain many features of the 
assembly, organisation and function of the non-centrosomal MT array in Drosophila oocytes 
and highlight the key role of diherential MT nucleation or anchoring at the cortex m 


RESULTS 

MTs in the oocyte are nucleated or anchored at the cortex psiiini and grow from the 
membrane into the volume. The cortical MT density follows a steep gradient along the 
posterior-lateral cortex from high density at the anterior corners to low densities at the 
posterior pole (Fig. ID, Materials and methods M2.1). In our model, we emulated these 
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features by selecting seeding points for MTs at random positions along the oocyte cortex, 
comprised of two parabolic, rotationally symmetric caps that capture the typical shape of 
a stage 9 Drosophila oocyte (Fig. lA, Materials and methods Ml). The density of seeding 
points decreases steeply along the posterior-lateral cortex to zero at the posterior pole, and 
decreases weakly along the anterior cortex towards the anterior centre (Fig. lA-B, Materials 
and methods M2.1). Each seeding point nucleates a MT polymer that grows until it either 
hits a boundary or reaches a target length imposed by the aging of MTs before catastrophe 
[2H] (Materials and methods M2.3). In total, we computed more than 55000 MTs for each 
realisation of a 3D wild-type cytoskeleton. 

A cross section through the computed MT meshwork (Fig. 1A,B) bears a striking resem¬ 
blance to confocal images of Tau- |29l [30] and EB-1- [16] tagged MTs, correctly showing the 
pronounced anterior-posterior gradient of MT density. Taking into account the orientations 
of all MT segments in the 3D volume, it also reproduces the experimentally measured direc¬ 
tional bias [IE] with 8.5% more MT segments pointing posteriorly than anteriorly (Fig. lA, 
inset). The directional bias in a 2D slice varies depending on the depth of the slice in the 
oocyte, ranging from 50.7% of posteriorly oriented MT segments at the cortex to 60.6% in 
the mid plane (Fig. IB, inset). This demonstrates that measurements in 2D confocal slices 
poorly characterise the fully extended 3D system. 

Central to understanding mRNA localisation is the question of MT orientations. Cal¬ 
culating the local vectorial sum of MT segments on a coarse-grained grid gives the local 
MT orientations in the computed cytoskeleton. Local orientations of an individual realisa¬ 
tion of the cytoskeleton (Fig. 1A,B) show poor global network order (Fig. 1C). However, 
cargo localisation in the oocyte does not involve only a single cytoskeletal realisation. MTs 
disappear within minutes in the presence of the MT-depolymerising drug colchicine that 
blocks MT growth by sequestering free tubulin dimers, indicating that the whole network 
turns over rapidly [151 [31] (V. Trovisco, personal communication). Thus, the oocyte sam¬ 
ples many tens to a hundred of independent MT organisations over the 6-9 hours of stage 
9. Summation of local orientations over an ensemble of 50 independent realisations of the 
computed MT network reveals a striking spatial partitioning of the cytoskeleton into several 
subcompartments (Fig. IF). At the anterior, the mean orientation points posteriorly, while 
MTs at the lateral sides on average point inwards toward the AP-axis. Counting the number 
of MT segments that contribute to each grid box in the ensemble also shows the distribution 
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of MT density (Fig. IH). In agreement with experiments (Fig. ID), the MT density exhibits 
a prononnced anterior-posterior gradient with highest valnes in the anterior corners, even 
when MT seeding is nniform on the anterior snrface. 

In the limit of a large ensemble, the MT polymers sample all possible initial directions 
and possible lengths at every point along the cortex. In this limit, we can test the predicted 
compartmentalised MT organisation by constrncting a second model in which MTs are 
represented as straight rods. The net orientation of the cytoskeleton at a given point inside 
the oocyte is then compnted as weighted snm of MT contribntions from each point along 
the bonndary (Materials and methods M3.1). 

With a shorter MT target length to compensate for the effectively longer length of straight 
rods compared to cnrved polymers, compntation of net MT orientations in the rod model 
shows a topology (Fig. IG) that conhrms the mean topology in the polymer model (Fig. 
IF). The three compartments of a posterior-pointing anterior section and two inwards- 
pointing lateral sections are effectively bonnded by separatrices (Fig. IG, orange arrows). 
Gargo molecnles that are transported on MTs to their pins ends move along the arrows 
and converge at the separatrices, eventnally leading to the posterior pole as the sole point of 
attraction in the entire oocyte volnme (the attractor). By contrast, cargo that is transported 
to MT minns ends moves opposite to the arrows and diverges away from the separatrices. 

This mean topology of the MT cytoskeleton is insensitive to the exact choice of the 
MT nncleation probability and to the choice of the MT length distribntion (Materials and 
methods M2.1, M2.3). It also remains nnchanged in a differently shaped oocyte geome¬ 
try in both the polymer model and the rod model (Fignre 1 - hgnre snpplement 1). In 
snmmary, both models show that even disordered non-centrosomal MT arrays can featnre 
well-dehned mean organisations, and that the MT cytoskeleton in oocytes is organised in a 
compartmental fashion. 

The snitably scaled local vectorial snm of MT segments (Fig. IG) for an individnal 
realisation of the polymer model (Fig. lA-B) is a vector held which represents active 
Kinesin-driven transport on the cytoskeleton (Materials and methods M2.5). In vivo dnring 
stage 9, the oocyte cytoplasm nndergoes slow cytoplasmic hows that are abolished in kinesin 
heavy chain mntants. This indicates that hows are driven by kinesin-dependent transport 
of an nnknown cargo throngh the viscons cytoplasm [26l |32] , thus making cytoplasmic hows 
a secondary read-out of cytoskeletal organisation [33]. Therefore, we next tested if our 
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computed polymer MT cytoskeleton can produce flows that are consistent with observed 
cytoplasmic streaming. 

Measurements of speeds of autofluorescent yolk granules in live stage 9 oocytes showed 
that mean flow velocities are slow (Fig. 2F). The physics of slow incompressible fluid flows 
u driven by forces / is described by the Stokes equations 

0 = —Vp + /iV^w + / , V • = 0 . (1) 

We make the simplest possible assumption that the forces / are proportional to the motor- 
velocity held, and use experimentally measured flow speeds to calibrate the scalar factor of 
proportionality. By solving the Stokes equations, we then computed the full 3D fluid flow 
field (Fig. 2A) corresponding to an individual realisation of the MT cytoskeleton (Fig. lA), 
and 2D cross sections through the 3D held (Fig. 2B) were compared to in vivo flow patterns 
visualised by particle image velocimetry (PIV, Fig. 2C). 

Despite large variability, both computed and in vivo flows show very similar patterns, 
generally being strongest in the anterior half of the oocyte, and weaker in the posterior half 
(Fig. 2B-C). Computed flows occasionally reach further into the posterior half than typically 
seen in PIV flow fields, thereby slightly overestimating the range of flows. However, except in 
rare cases, computed flows do not reach the posterior pole. This result is largely independent 
of the presence of the oocyte nucleus which occupies less than 2.5% of the oocyte volume, 
even for small oocytes at the beginning of stage 9 (Figure 2 - figure supplement 1, Materials 
and methods M5.2, M5.3). Thus, computed flows appear consistent with observed slow 
cytoplasmic streaming. 

We tested next if our computed cytoskeleton in combination with the derived cytoplasmic 
flows and diffusion can account for dynamic mRNA transport and localisation. We described 
the mRNA distributions as continuous concentration fields, oskar mRNA is assumed to 
reside in either one of two states: the Kinesin-bound state with concentration c;,, in which 
cargo is transported actively on the cytoskeleton-derived motor-velocity held Vm (Fig. 1C), 
or the unbound state with concentration Cu, in which cargo is transported by the cytoplasmic 
flows u (Fig. 2A) and diffuses with diffusion constant D. Cargo can exchange between both 
states by chemical reactions, thereby resulting in the reaction-advection-diffusion equations: 
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Parameter values were constrained with experimentally measured values (Materials and 
methods M4.2). Furthermore, throughout the simulation we cycled through the pairs of 
fluid flow u and motor-velocity helds (Fig. 3) to account for the dynamic nature of the 
MT network [15] and flow patterns, which change over time scales of minutes compared to 
the 6-9 hours during which oskar mRNA is localised (Materials and methods M5.1). 

Starting from an initial diffuse cloud of oskar mRNA in the centre of the oocyte (Fig. 
3A,D,E), simulations show that the mRNA quickly concentrates in the centre before forming 
a channel from the centre of the oocyte towards the posterior pole (Fig. 3E). Formation of 
this channel reflects the locally inwards orientation of MTs in the posterior half of the oocyte 
(Fig. 1E,F). This transient state closely resembles transient patterns of fluorescently-tagged 
oskar mRNA during the transition between stages 8 and 9 (Fig. 3B), thereby supporting 
the notion that our computed MT cytoskeleton correctly captures key aspects of the in vivo 
MT network. 

Upon reaching the posterior cortex, oskar mRNA is translated to produce Long Osk [3l] , 
which is involved in anchoring oskar mRNA. However, oskar mRNA localises normally at 
stage 9 when anchoring is disrupted [3ll|35]. Anchoring is therefore not necessary at stage 9, 
and we did not include it in our model at this point. Despite the lack of any anchoring, oskar 
mRNA forms a posterior crescent after 1.5 hours and becomes highly concentrated by the 
end of the simulation (Fig. 3F), correctly capturing observations of in vivo localisation (Fig. 
3C). A 4/im thick slice at the posterior pole contains 12.5% of total cargo in the oocyte (Fig. 
3F), showing that continual transport combined with slow diffusion is sufficient to reach and 
maintain high concentrations of mRNA. 

In our description of transport (Eq. (2)), cargo acts as an unspecihc passive tracer. 
Passive tracers merely follow the transport helds and u, and these two helds must contain 
all information about the localisation sites. We therefore asked which held contains most 
information about localisation. We hrst tested localisation in the absence of cytoplasmic 
hows by setting the huid how held identical to zero u = 0. This situation is similar to but 
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more severe than in slow Kinesin mntants [26]. Starting from the central oskar mRNA clond 
(Fig. 3D), we fonnd that oskar mRNA still localises to and forms a concentrated crescent at 
the posterior pole (Fig. 3G). A 4/im thick slice off the posterior pole contains 11.3% of total 
available oskar mRNA cargo, only marginally less than localisation with cytoplasmic flows 
present. This snggests that the cytoskeletal transport alone localises the majority of mRNA. 
Cargo localisation in the absence of cytoskeletal transport can be tested by either setting the 
motor-velocity held identical to zero Vm = 0, or by setting the binding constant ki, to zero. 
In either case, cargo disperses thronghont the oocyte withont forming a central channel and 
completely fails to produce any posterior crescent. We further asked whether the how could 
localise cargo to wild-type levels if an anchor captured cargo at the posterior. To this end, 
we replaced the bound state in the model by an anchored state to which cargo can bind only 
at the very posterior boundary. Binding to the anchor occurs with a reaction rate constant 
that is ten times larger than the reaction rate constant kh used for oskar mRNA binding to 
the MT cytoskeleton (Materials and methods M4.2). Moreover, anchored cargo cannot be 
released, thus making this anchor a perfect sink with stronger trapping properties than any 
realistic anchor. Under these idealised conditions, oskar mRNA accumulates slightly at the 
anchor (Fig. 3K). However, a 4/im thick posterior slice contains only 2.7% of total cargo 
in the oocyte, about 78% less cargo than in the wild-type simulations, and only 25% more 
cargo compared to the purely diffusive case, which localises 2.2% of cargo. Thus, even under 
the most favourable conditions, the cytoplasmic flows cannot localise mRNA to wild-type 
levels, arguing against a mixing and entrapment mechanism for oskar mRNA localisation 
at stage 9. Despite the small fraction of actively transported cargo (13%), motor-driven 
transport is both necessary and sufficient to account for oskar mRNA localisation. 

In contrast to oskar mRNA, bicoid mRNA is believed to be transported by Dynein pra]- 
EU, but other parameters such as the fraction of bound cargo and the active transport speeds 
are similar to oskar mRNA [36] (V. Trovisco, personal communication). To test whether 
the proposed cytoskeletal organisation also captures transport and localisation of injected 
bicoid mRNA, we inverted the directions of the motor-velocity helds (Fig. 1C) to account 
for the minus end directed transport by Dynein instead of the plus end directed transport 
by Kinesin. Other parameters including cytoplasmic flow helds remained unchanged. 

Closely matching experiments with injected naive bicoid mRNA im. simulations show 
accumulation of bicoid mRNA at both posterior-lateral sides of the oocyte when placed 


initially in the posterior half (Fig. 3H, inset), or accumulation at the anterior and ventral 
cortex when initially placed in the anterior-ventral area (Fig. 31, inset). Localisation of 
bicoid mRNA is also virtually identical in the absence of cytoplasmic streaming. Thus, 
splitting of a bulk amount of injected bicoid mRNA occurs when the RNA is placed on the 
border (separatrices) between two diverging subcompartments of the MT cytoskeleton, each 
one transporting part of the cloud towards the adjacent cortex. The agreement between 
simulations and experiments therefore further supports an on average compartmentalised 
MT cytoskeleton, and naive bicoid mRNA, like oskar mRNA, unspecihcally traces out the 
MT cytoskeleton. 

In contrast to naive bicoid mRNA, endogenous bicoid mRNA is transported into the 
oocyte after being transcribed in neighbouring nurse cells where it also becomes modihed 
in presence of the protein Exuperantia. Endogenous bicoid mRNA can be approximated by 
so-called conditioned bicoid mRNA that is hrst injected into the nurse cells, then sucked out 
and injected into the oocyte. Conditioned bicoid mRNA localises specihcally to the anterior 
surface within 30 minutes, even if this surface is not closest to the injection site HU. In the 
model, starting from an initial distribution at the anterior surface where endogenous bicoid 
mRNA enters the oocyte, cargo quickly moves to the anterior cortex before concentrating 
in the anterior corners after several hours of simulated time (Fig. 3J). This resembles 
the anterior ring-like localisation of bicoid mRNA in wild-type stage 9 oocytes EBEZ!, and 
occurs independently of cytoplasmic flows. However, the model does not reproduce transport 
specihcally to the anterior surface when injection is further away from the anterior. This 
suggests that Exuperantia activity somehow either enables bicoid mRNA to move along an 
unidentihed population of MTs naE] that are not included in the computed cytoskeleton, 
or that it allows localisation of the mRNA by an unknown MT-independent mechanism. 

We have shown so far that cortical gradients of MT nucleation are sufficient to assemble a 
functional, compartmentalised MT cytoskeleton that successfully localises oskar mRNA and 
naive bicoid mRNA. Next, we asked if cortical MT nucleation can also produce cytoskeletal 
organisations that explain mutants with partial or complete polarity defects. Mutations 
interfering with the posterior follicle cells that surround the oocyte, for example, can disrupt 
the proper positioning of mRNAs [381 1^ . Specihcally, follicle cell clones (ras^‘^^°'’) that are 
adjacent to only one side of the oocyte posterior repel cargo from that side of the cortex 
(termed clone adjacent mislocalisation jlQ]), likely due to an altered MT organisation. To 
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mimic this phenotype, we used an ensemble of wild-type cytoskeletons and artihcially added 
MT nucleation sites to a patch on one side of the posterior pole. Simulations of cargo 
transport with diffusion, motor transport and flows show that oskar mRNA is repelled from 
this patch and localises to the adjacent posterior boundary, in agreement with experiments 
(Figure 4 - hgure supplement 2). 

Mislocalisation of mRNAs also occurs in mutants of the polarity protein PAR-1, which 
acts to suppress MT nucleation at the posterior pole of the oocyte |H1112]. In strong par-1 
hypomorphs with strongly reduced PAR-1 expression, MT minus ends occupy the whole 
cortex including the posterior pole [I6], thereby causing oskar mRNA to mislocalise to a 
dot in the centre of almost 90% of mutant oocytes [22] . To test this phenotype in the model 
we distributed MT seeding points uniformly on the posterior surface (Fig. 4G, inset) while 
keeping the seeding density on the anterior surface constant. The anterior-posterior gradient 
of MT density then vanishes and the directional bias of MT segments evens out (Fig. 4G, 
ensemble-averaged 3D-bias: 50.1%:49.9%). The ensemble-averaged local orientations of MTs 
show that MTs point towards a single focal point in the centre of the oocyte (Fig. 4G) that 
acts as a stable hxed point of the system. Gomputation of 3D flow helds for each realisation 
in the new ensemble (3D mean: 13.5 nm/s; 2D mean: 14.5 nm/s, N=50) and simulations of 
transport with and without cytoplasmic flows show oskar mRNA concentrating in a cloud in 
the centre of the oocyte (Fig. 4G). oskar mRNA accumulation to a central dot also remains 
unchanged if the posterior MT seeding density is decreased slightly (Fig. 4H, inset) as might 
be expected for hypomorphs that do not abolish PAR-1 expression completely. 

Interestingly, bicoid mRNA simulated either with slightly decreased (Fig. 4H) or with 
uniform (Fig. 4G) posterior MT nucleation not only localises to the anterior corners as in 
wild-type but also enriches at the posterior pole (Fig. 4D, arrowheads), despite injection 
close to the anterior (Fig. 4D, inset). This matches experiments showing bicoid mRNA 
mislocalisation to both anterior and posterior in gurken, torpedo and cornichon mutants 
(EGFR mutants) in which posterior follicle cells fail to signal a MT reorganisation in the 
oocyte [38l l39] . 

Failure of the external signal to the oocyte is thought to have the same effect as par-1 
hypomorphs, but previous experiments with bicoid mRNA in par-1 hypomorphs resulted 
in ambiguous localisations [111113]. Here, using confocal fluorescence in-situ hybridisation 
(FISH), our experiments show that bicoid mRNA localises primarily to the anterior corners 
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and to the posterior pole in strong par-1 hypomorphs (Fig. 4L), thereby mirroring EGFR 
mntants and agreeing with the model predictions. Thns, (near-)nniform MT nncleation along 
the posterior and lateral cortex is snfficient to explain the polarity phenotypes observed for 
oskar and bicoid mRNAs in gurken, torpedo and cornichon mntants and in strong par-1 
hypomorphs. 

In 10% of strong par-1 hypomorphs and 30% of weak par-1 hypomorphs with weakly 
rednced PAR-1 expression, oskar mRNA is only partially mislocalised, combining a wild- 
type-like posterior crescent with a centrally mislocalised dot [22]. We therefore modnlated 
the MT nncleation along the posterior-lateral cortex to test if this is snfficient to prodnce 
intermediate cytoskeletal organisations between wild-type (Fig. 4F) and strong par-1 hypo¬ 
morphs (Fig. 4H). 

Increasing MT nncleation laterally towards the posterior pole (Fig. 4G, inset) again 
creates a stable focns in the centre of the oocyte. However, its domain of attraction does 
not cover the entire oocyte. Instead, a small basin of attraction towards the posterior pole 
persists (Fig. 4G, red), separated from the basin of attraction of the stable focns by an 
nnstable saddle-node point. Fnrther increasing the seeding density towards the posterior 
pole shows that the pair of stable and nnstable hxed points move fnrther apart (Fignre 4 - 
hgnre snpplement 1) nntil the nnstable point no longer resides inside the oocyte geometry, 
thereby giving rise to the strong par-1 hypomorph topology (Fig. 4H, I, J). Gargo simnlated 
on cytoskeletons with stable and nnstable hxed points and their corresponding how helds can 
split into two separate accumnlations, hrst at the posterior pole and second in a dot along 
the AP-axis (Fig. 4B). This agreement with experimental resnlts snggests that intermediate 
levels of PAR-1 at the posterior lead to a cytoskeleton with a potential barrier between the 
oocyte centre and its posterior pole. It also shows that modnlation of MT nncleation along 
the posterior-lateral cortex alone is snfficient to captnre this. Becanse the potential barrier 
only ahects the bonnd state, dihnsion and hows can aid posterior localisation in this context 
by pnshing nnbonnd cargo across the nnstable point and increasing the amonnt of oskar 
mRNA that reaches the posterior pole. In snmmary, variations of MT nncleation along 
the cortex not only explain mRNA localisations in wild-type bnt can also acconnt for the 
mislocalisations of bicoid and oskar mRNAs in polarity mntants. 

The creation of stable and unstable hxed points in the transition from wild-type to 
strong par-1 hypomorph topology constitutes a classical saddle-node bifurcation (Fig. 4K- 
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M, Materials and methods M3.3) with the cortical MT seeding density as the bifnrcation 
parameter. Interestingly, in both the polymer model and the rod model, the mean length 
of MTs acts as a second bifnrcation parameter. For example, for snfficiently short MT 
hlaments (Fig. 4K), few MTs from the anterior can reach and contribnte to orientations 
in the posterior half of the oocyte. In the posterior half, MTs from the lateral side either 
point towards the anterior, thereby combining with anterior MTs to create a stable node, 
or towards the posterior, to form a domain of attraction at the posterior pole. Therefore, 
for snfficiently weak contribntions from anterior MTs, MTs from the lateral sides create the 
unstable tipping-point. 


To understand more generally how the MT lengths and the posterior-lateral distribution 
of seeding points together influence the three different cytoskeletal topologies (Fig. 4E- 
G) we computed the parameter space of the straight rod model (Fig. 5). If MTs are 
absent at the posterior pole, the wild-type topology (Fig. 5D i, vi) covers a large fraction 
of parameter space (Fig. 5A, red) for sufficiently long MTs (large e) and cortical MT 
nucleation gradient (large kp). The creation of two fixed points splits the oocyte into 
different domains of attraction (Fig. 5A, blue, D i-ii) for either shorter MTs (Fig. 5A, 
vertical dashed arrow; Figure 4 - figure supplement 1 I,N) or for increasingly uniform MT 
nucleation (Fig. 5A, horizontal dashed arrow; Figure 4 - hgure supplement 11,H). For almost 
uniform nucleation along the entire posterior-lateral cortex, the unstable hxed point exits 
the geometry, leaving behind only a single stable focus (Fig. 5A, green, arrowhead, D iii-iv). 
If some MT nucleation is allowed at the posterior pole this region of strong par-1 hypomorph 
topology expands substantially (Fig. 5B, C, green). Interestingly, if MTs are sufficiently 
long, MTs from the lateral and anterior cortex can reach and overpower those from the 
posterior pole, thus restoring wild-type-like cytoskeletal topology that allows oskar mRNA 
localisation. Therefore, a low level of posterior nucleation of MTs does not necessarily lead 
to mislocalisation of oskar mRNA. Instead, and in addition to the distribution of cortical 
MT nucleation, the length of MT filaments emerges as another key regulator for polarisation 
and function of the non-centrosomal network. 
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DISCUSSION 


Non-centrosomal MT networks represent a large, yet poorly understood class of MT ar¬ 
rangements that often fulfil specialised functions [H]. Non-centrosomal MTs are frequently 
aligned in parallel, thereby forming linear arrays as in epithelia or neurons [H El S] • The 
non-centrosomal MT cytoskeleton in Drosophila oocytes was also hrst hypothesised to form 
a highly polarised linear array with MTs growing from the anterior surface towards the 
posterior pole [iHl EH] . Instead, the MT cytoskeleton is a complex network without a vis¬ 
ibly discernible polarity along the anterior-posterior axis [la na Eg, and its organisation 
remained ambiguous 

We showed that in two different theoretical models, in which MTs grow from the oocyte 
cortex into the volume, the cytoskeleton features an on average compartmental organisation 
and is therefore more ordered than previously thought. In combination with cytoplasmic 
flows and diffusion, this cytoskeletal organisation successfully reproduces localisations of os- 
kar and naive bicoid mRNAs in wild-type oocytes. Such a spatially varying, compartmental 
organisation suggests that a single statistical measure of polarity, derived by averaging data 
from the entire oocyte mils], may not be a sufficient metric to characterise mRNA locali¬ 
sation. 

Cytoplasmic flows may generally contribute to mRNA transport. Flows at later stages of 
oogenesis are fast and well-ordered, and occur concurrently with a late-phase enhancement 
of oskar mRNA accumulation m- This is consistent with a contribution of flows to mRNA 
localisation [211 I2S]- At stage 9 of oocyte development, which is the focus of the work 
presented here, flows are slow and chaotic. At this stage, mutants with reduced flows 
showed oskar mRNA localising normally [26], thus suggesting that slow and chaotic flows 
may not play the same role as in late oogenesis. However, flows in these mutants were 
merely reduced rather than abolished completely, and the measurements underestimated 
flow speeds [21], hence leaving the interpretation of flows unclear. We here hnd that the 
effects of slow cytoplasmic flows on mRNA transport are negligible, and that cytoskeletal 
transport alone is sufficient for localisations of oskar and naive bicoid mRNAs. In this view, 
slow cytoplasmic flows arise primarily as inevitable physical byproduct of active motor-driven 
transport on the cytoskeleton rather than as an evolutionarily selected trait. This appears 
to mirror hndings in the C. elegans zygote in which P-granules segregate by dissolution and 
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condensation rather than via transport by cytoplasmic flows [H] . 

Interestingly, neither MT-based transport alone nor combined with cytoplasmic flows and 
diffusion is sufficient to reproduce the anterior localisation of nurse cell conditioned bicoid 
mRNA irrespective of the position of injection into the oocyte im. This Exuperantia- 
dependent mechanism may rely either on an unobserved population of MTs in the oocyte 
that can be specifically recognised by conditioned bicoid mRNA [min], or on an unknown 
MT-independent mechanism. 

The central finding of our work is that gradients of cortical MT nucleation are sufficient 
for the assembly of a functional compartmentalised MT cytoskeleton in wild-type oocytes. 
While many non-centrosomal MT arrays are linear and emphasise questions about estab¬ 
lishment and maintenance of parallel filament orientations l50] , our result stresses the 
need to understand gradients in MT nucleation as an alternative strategy for the assembly 
of functional non-centrosomal arrays. Whether MTs in Drosophila oocytes are differentially 
nucleated along the cortex itself or created elsewhere and then differentially anchored at the 
cortex remains an interesting open question, but both scenarios are compatible with our 
model. Understanding how this gradient is established will therefore depend on discovering 
how PAR-1 regulates MT interactions with the cortex. 

Cortical MT gradients cannot only account for the wild-type cytoskeletal configuration 
but also for the phenotypes observed in a hierarchy of par-1 hypomorphs. Modulation of 
MT gradients along the posterior-lateral cortex alone are sufficient to explain the splitting 
of oskar mRNA between the centre and the posterior pole of the oocyte [22] via a saddle- 
node bifurcation, suggesting that the anterior and posterior-lateral surfaces of the oocyte 
are functionally decoupled. Generally, bifurcations in temporal behaviour govern important 
qualitative transitions in many biological systems, such as the lactose network in E. coli 
[SI]) cell cycle in yeast [S2]) ^md collapses of bacterial populations [SB]- In Drosophila, one 
important qualitative change is the temporal transition between stages 7/8 and 9. During 
this transition the MT cytoskeleton reorganises from uniform nucleation around the cor¬ 
tex and oskar mRNA in the centre to the anterior-posterior MT nucleation gradient with 
oskar mRNA at the posterior. It is therefore tempting to speculate that the sequence of 
PAR-1 mutants and the underlying bifurcation represent static snapshots of this dynamic 
developmental transition in wild-type. 

In conclusion, the present work provides a model that describes the assembly, organisation 
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and function of the non-centrosomal MT array in Drosophila oocytes and directs future 
attention to the molecular mechanisms that enable differential MT nucleation or anchoring 
at the cortex. 

MATERIAL AND METHODS 

Ml Coordinates and oocyte geometries 

For calculation of the microtubule (MT) cytoskeleton, the anterior-posterior (AP) axis 
is aligned with the ; 2 -axis of a cartesian coordinate system, and results are later shifted and 
rotated to align with the x-axis for visualization, for subsequent computations of cytoplas¬ 
mic flows and simulations of cargo transport. Dimensional coordinates cover the ranges 
X G [—L, L], y E [—L, L] and G [0, Zq L] with length scale L = 50 ym. After nondimension- 
alization with scale L, coordinates are ranged as x' G [—1,1], y' G [—1,1] and z' G [0,2:o]. 
Similarly, nondimensionalization of the shape parameter kp in the MT seeding density leads 
to k'p = kp/L, and all primes will be dropped subsequently. Methods hgures are shown 
in nondimensional spatial coordinates, and reported parameter values are nondimensional 
unless noted otherwise. 

The 3D geometry of a typical stage 9 Drosophila oocyte is dehned as two parabolic, 
rotationally-symmetric caps (Fig. 1). Anterior [i = A) and posterior [i = P) parabolic caps 
are parameterized in cylindrical coordinates p = G [0,1] and 0 G [0, 2ti) as 

= pep + zl{l- (3) 

where Cp = (cos(0), sin(0), 0) and = (0, 0,1) are the unit vectors in radial and 2 ;-direction, 
respectively. The line element along the parabola is calculated as 

dsj = ^ dp = \Jl + ( 24 p)^dp 

with arclength 

s*(p) = ^ \/l + (24 ^ sinh-i (24 p) (4) 

dehned such that Si(p = 0) = 0 denotes the tip of the parabolic cap while Si{p = 1) = Sg 
denotes the distance from the tip to the anterior corners (Fig. 6 A,G), and hence 0 < <Sj(p) < 
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Sq. The surface element for parabolic caps (in units of L) is given by 




dcTi dcTi 

dp ^ d(j) 


P\f^ 


( 24 p)^dpd 0 , 


(5) 


with total surface area 


S, = 


TT 


6 {z^y 


([1 + ( 24 )^]=‘-'^-i) 


The inward pointing normals for posterior and anterior caps are given by 

np{p, 0) = 

nA{pA) = 


2z^ pBp + e^ 

a/i + yiz^py 

2zq pCp + 

Vi + i:^4p? ■ 


( 6 ) 

(7) 


A special case arises for Zq = Q when the anterior parabolic cap becomes a flat disc. In this 
case, the arclength (eq. (|^) reduces to sa{p) = p, and the surface area (eq. (|^) simplihes 
to Tja = tt. The inwards pointing normals for this case are given by np (eq. 0 ) for the 
posterior cap and for the anterior disc, respectively. 

To investigate the robustness with respect to changes in oocyte shape, we test our model 
for the cytoskeleton, cytoplasmic flows and mRNA transport in two different geometries for 
a stage 9 Drosophila oocyte. Geometry-1 is used as standard geometry, and is comprised 
of two parabolic caps (eq. (|^) with Zq = 1.48 for the posterior cap and z^ = 0.2 for the 
anterior cap (see Fig. 1). This results in a length of the anterior-posterior axis (AP-axis) of 
Zq — z^ = 1.28 with an aspect ratio of 1.56 that qualitatively resembles a typical stage 9 
oocyte. Geometry-2 is tested as alternative geometry, and consists of a posterior parabolic 
cap with Zq = 1, and an anterior flat disc with z^ = 0 (Fig. 1 - hgure supplement 1), giving 
an AP-axis length of 1 and aspect ratio of 2. We hnd that results are robust with respect 
to this change in geometry. 


M2 Polymer model for MT cytoskeleton 

M2.1 MT nucleation probability 

The hrst step in the computation of the MT cytoskeleton is the generation of MT seeding 
points that are randomly positioned along the oocyte membrane. We dehne the following 


16 














probability density along the arclength Sj(p) of the anterior and posterior parabolic caps 



(8) 


with normalization over the oocyte snrface 



<l'LAPT.{sA{p)\hQ,kA) 


dSpps(sp(p)|h^, fcp) , 


(9) 


where the snrface elements were specihed in eq. (|^ and the arclength was given in eq. (|^. 
For the fnnctional form of ps("S(p)|/io, k) we use the expression 


Ps(s(p)|ho,fc) = ho + (1 - ho) 1 + 




( 10 ) 


+ s{pY ’ 


wherein the index i was dropped for simplicity. Note that for the maximum arclength s = sq 
it holds that pe(s = Sojho, k) = 1. 

ho G [0,1] and k G (0, cx)) represent two parameters governing the shape of the probability 
density. For small values h so the nucleation density approaches a Hill function with Hill 
coefficient n = 2 and depth ho 



( 11 ) 


In the opposite limit h 3> so, the probability density approaches a parabola with depth ho 
(Fig. 6B,H, red curve) 



( 12 ) 


thus creating a pronounced gradient of MT nucleation from the pole (s = 0) to the corners 
(s = So). 

MT nucleation along the arclength can increase from the gradient to a homogeneous 
nucleation along the arclength in two different ways: (i) either by increasing nucleation from 
the corners laterally towards the posterior pole, or (ii) by increasing nucleation uniformly at 
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the posterior pole. Each shape parameter ho and k accounts for one of these possibilities. 
Increasing the parameter ho leads to shallower gradients by increasing the probability for 
nucleation at the pole from 0 for ho = 0 to the value for homogeneous nucleation for ho = 1 
(Fig. 6H). Contrary, decreasing the parameter k leads to narrower regions of low nucleation 
probability by decreasing its width laterally (Fig. 6B). 

To compare to experimental data, we stained a-tubulin in hxed stage 9 wild-type 
Drosophila oocytes and extracted fluorescence intensity prohles along the cortex from the 
posterior pole to the anterior corners (Fig. 7A,B). An ensemble of smoothed fluorescence 
prohles shows steep increases in cortical MT density from posterior to anterior, albeit with 
high variability. Still, a quadratic ht of the ensemble shows that a parabolic function is 
a plausible representation of the increasing MT density (Fig. 7C). Several oocytes show 
huorescent prohles that exhibit a local maximum at the posterior. This likely stems from 
out of plane huorescence due to the geometric shape of the oocyte pole, and we neglect 
this minor ehect here. At the anterior surface, MT density tends to be more homogeneous 
with only shallow increases in cortical MT density towards the corners. In summary, for the 
posterior cap we choose a probability density with parabolic shape that reaches zero at the 
posterior pole = 0,kp = 20 (Fig. 7D, bottom), whereas for the anterior cap we chose a 
density with parabolic shape but only shallow depth = 0.8, kp = 20 (Fig. 7D, top). 


M2.2 Generation of seeding points 

Discrete seeding points that are randomly positioned along the oocyte membrane accord¬ 
ing to the probability density Ps(-Si(p)|ho, ki) (eq. ([^) can be achieved by inverse transform 
sampling. To this end, the cumulative distribution function on the anterior {i = A) or 
posterior {i = P) cap can be computed as 

p2TT pp' 

C{p'\hl,ki) = / / dSips(si(p)|ho,fci) . (13) 

Jo Jo 

Note that C{p' = l\hQ,ki) ^ 1 due to normalization across the total oocyte surface rather 
than across each individual cap. Therefore, random seeding points on each cap can be 
drawn by applying the inverted cumulative distribution function to a random sample that 
is uniformly distributed over the interval [0,C(p' = l|ho,/ci)]. In this formulation, the total 
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number of seeding points is iV = Na + Np (see Fig. 6). The ratio of anterior (Na) and 
posterior (Np) seeding points to be drawn must equal the ratio of total probability for a 
point to fall on the anterior and posterior caps, respectively. For a hxed value Na, this 
results in 


^ JdEppj:{sp{p)\h^,kp) 

^ ^ J dJ:APj:{sA{p)\ho,kA) ' 


(14) 


Here, we generated seeding points in an equivalent way by renormalizing the probability 
density (eq. (|^) on each cap individually as 


ki) 


1 

A 


Pj:{si{p)\hl, ki) , 


(15) 


wherein the new normalization factor A is obtained by integration over only one parabolic 
cap 

A{hl, ki) = J dSjps(si(p)|ho, ki) . (16) 

Hence, the corresponding cumulative distribution function 


cip'Kh) 



dSiF’s(si(p)|ho, ki 


(17) 


covers the range [0,1] and seeding points are obtained as C ^{u) for a random sample u with 
uniform distribution u G [0,1]. 


For parabolic caps, the normalization (eq. imi) can not be calculated analytically. There¬ 
fore, the probability density (eq. (Hi) and its cumulative distribution function (eq. ([I7|) 
were evaluated and inverted numerically. For the special case of a flat anterior {i = A) disc 
in the alternative geometry-2, the normalization can be solved in closed form as 




A{hQ,kA) = J d(p J dp ppY.{p,hQ,kA) = (18) 

TT [(1 - h^)k\ + 1 - (1 - h^)k\{k\ + 1) ln(l + l/k\)] . 


The corresponding cumulative distributions function evaluates to 


/ 27 r np' 

d(j) j dppV^{p,hQ,kA) = 

p^ [(1 - K)k\ + !]-(!- K)k\{k\ + 1) ln[l + (p/fc^)^] 
(1 - h^)kl + 1 - (1 - K)kl{kl + 1) ln[l + [l/kAf] 


(19) 
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which was numerically inverted to generate randomly positioned seeding points on the an¬ 
terior disc. 


Due to the new individual normalization, the values of anterior and posterior probability 
densities differ at the contact line between anterior and posterior caps. Consider two rings 
with areas dAj along the contact line (s = Sq) on the anterior {i = A) and posterior {i = 
P) cap which contain dNi = NiVj:{s = SqI/iq, fcj) dA* seeding points. Setting the number 
of anterior points fixed, and enforcing the point densities to be equal in the corner rings 
dA^^/dA^ = dNp/dAp for equal ring areas dA^ = dAp allows to compute the number of 
posterior seeding points as 


Np = Na 


= Na 


V^{s = s^\h^,kA) 
= s^\h^,kp) 
Aih^^kp) 
Aih^kA) ’ 


( 20 ) 


where the second line follows because ps^Si = SQ|h.Q,/cj) = 1. Inserting the expression for 


A{hQ, ki) (eq. (16)) shows that forcing seeding point densities on the anterior and posterior 
caps to be equal at the anterior corners (eq. (|20|) is identical to chosing points according to 
the ratio of total probabilities on the anterior and posterior caps (eq. (14)). Thus, chosing 
points from ps normalized over the entire surface (eq. or from Vp, normalized over each 
cap (eq. 15) is equivalent. 


Alternative seeding densities Our model is not sensitive to the exact choice of the proba¬ 
bility density for MT seeding points, and here we define alternative seeding densities that are 
mathematically tractable. Specifically, as randomly positioned seeding points are generated 
by operating the inverse CDF on a uniformly distributed sample u G [0,1] and (j) G [0, 27r), it 
is convenient to chose the seeding probability such that the CDF can be inverted analytically. 


We first note that writing out explicitly the CDF for a rotationally-symmetric probability 


density function p{si{p)) along the arclength (eq. (13)) leads to 


j^ 2 ^^ rp , -^- 

c(p') = <i<t> dpp(Si(p))pgi + (2zip)'- 


( 21 ) 


For the case of a gradient of MT seeding points along the arclength of a parabolic cap we 
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can define the density to be 


Pisiip)) =Pn{p\n,ZQ) 


{n + 1) p" ^ 

27r y/l + ( 24 p)' ' 


( 22 ) 


From eq. (21) this resnlts in a CDF u = Cn{p') 
as p' = C~^{u) = 


which can be analytically inverted 


For the case of a nniform distribntion of seeding points along the arclength, the probability 
density is constant and given by the inverse snrface area of the parabolic caps (eq. 


p(».(p)) =pM^o) 



(23) 


For Pu{p\zq) we hnd the CDF 


u = Cu{p') 


(1 + (HpTf -1 

(1 + (24)2)=''" - 1 ’ 


which can again be inverted analytically to give 


(24) 


c-‘(«) 



(1 + (24)2)=/" - 1 



(25) 


For the special case of a flat disc, the nniform distribntion of seeding points simplihes the 


CDF to 




dp-p 

TT 


P 


(26) 


with inverse p = C^^{u) = ^/u. 


For a wild-type cytoskeleton, we impose a nniform MT density Pu{p\zq) along the anterior 
cap and a gradient of MTs p„(p|n, Zq) along the posterior cap. At the contact line between 
anterior and posterior caps, the seeding densities need to be matched by adjnsting the total 
nnmber of posterior seeding points Np. A ring on the anterior cap at the anterior corners 
contains dAhi = {Na/^a) dA^ seeding points, while a ring on the posterior cap at the 
anterior corners contains dA'p = NpPn{p = l|n, Zq) dAp points. Using eq. (22) and eq. § 
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and forcing densities to be eqnal diV^/dA^ = dNp/dAp resnlts in 


^APn{p = 1 ) 

i(2z§YN^ sJl + (2z^Y 
" +1 (1 + {2z^yf''^ - 1 ' 


( 27 ) 


For the alternative geonietry-2 comprised of a posterior parabolic cap and flat anterior disc, 
matching of seeding densities resnlts in 


Np = 


2Na 
n + 1 


\/^ + . 


(28) 


Using these alternative seeding densities does not qnalitatively change the behavior of the 
model. 


M2.3 MT growth 


We describe MTs as persistent random walks, i.e. a chain of straight segments of constant 
length A that show some flexibility in their relative orientations. The orientation /ij of the 
f-th segment is drawn from a von Mises-Fisher probability distribntion on a 2D sphere 




K 

47rsinh(K) 


K. pji • fli — i 
^ 5 


(29) 


where k is the concentration parameter aronnd the orientation p,i-i of the previons segment 
(see sec. 2.5 for parameter valnes). The hrst orientation p,i is drawn from the nniform 
angnlar distribntion on a 2D sphere {k = 0) where orientations are rejected if pointing 
locally ontwards of the geometry fiA,p{p, (f>) ■ pLi < 0. 

Polymers stop growing if they either enconnter a bonndary, or if they reach a pre-dehned 
target length I drawn from a probability distribntion. The probability density is based on 
the experimental hnding that MTs in-vitro nndergo a three-step aging process leading to 
catastrophe, resnlting in catastrophe lengths (/c) that follow a Gamma distribntion [28] . 
From this catastrophe length distribntion, the length distribntion observed in an ensemble 
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of growing MTs without shrinkage was calculated as 


0r(%, A) = 


r(n,//A) 
n A r(n) 


(30) 


where r(n, /) is the incomplete Gamma function, r(n) is the Gamma distribution, n = 
3 is the number of steps to reach catastrophe and A is the step length (equivalent to a 
transition time per step for constant growth velocity, see [28]). We set the expectation value 
= 3,A)] = 2A = as fraction e of the anterior-posterior axis length, hence 


A = e (; 


/ 2 . 


Pseudo-random numbers distributed according to this probability density (eq. (30)) were 
generated by inverse transform sampling of the cumulative distribution function 


<J.r(i|n = 3. A) = 1- — -e-V-', 


(31) 


Similar to the MT seeding densities, the model is not sensitive to the exact choice of the 
MT length distribution. For example, using the one-parameter exponential distribution 

= (32) 

as an alternative length distribution and again regulating the expectation value E [0e(^|A)] = 
A = e (z(f — Zq) of MT lengths via a parameter e does not qualitatively change our results. 


M2.4 MT persistence length 

The stiffness of a polymer is commonly characterized by its persistence length P, which 
is dehned as the decay length of tangent-tangent correlations. To calculate the persistence 
length of MT random walks, we hrst dehne the end-to-end vector for a polymer with Ng 
segments of length A and given initial orientation fii as 

Ns 

i?(W|Ai) = Ai • (33) 

i=l 
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Orientations are drawn from a von Mises-Fisher distribntion (eq. (29)) with mean 


{x\p,) = j X K)dx = 


tanh(K) K 


fl = a p,. 


(34) 


In order to calculate the persistence length of these polymers, we hrst consider the mean 
orientations of the second and third MT segments 


(A2|Ai) = y p2f{P'2\fll,l^)dp2 = (^fll, 

{pslpl) = / As JJ/(/ii|/ij_i,K)d/ii = cr^Al , 

d i=2 

and we hnd in general the relation 


Ns 


/ Miv, Jf'/“i- 


(35) 


i=2 


Hence, using eq. (35) and shifting indices, the expectation value of the end-to-end vector 


{R{Ns\pi)) is given by 


Ns Ns-1 , _ TV, 

(AilAi) = A V] cr^Ai = A—- pi. (36) 

^^ ^^ i — CT 

i=l n=0 

For stiffness parameter k —)■ 0 (a —)■ 0), the von Mises-Fisher distribution becomes a uniform 
distribution on a 2D sphere. In this case, the polymers curl up with mean end-to-end distance 
of a single segment 

lim (R(Nsjpi)} = lim (R(Nsjpi)} = X pi- (37) 

K—^-O a—^0 

In the opposite limit k —)■ cx) (a —)■ 1), the polymers become stiff rods and are extended to 
their maximum possible length 


lim {R{Ns\pi)) = lim {R{Ns\pi)) = XN^ pi, (38) 

R^OO cr^l 

thereby showing that the polymers interpolate between an undirected random walk and 
completely straight rods. 

The persistence length P is dehned as the spatial decay length of the correlations of 
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polymer segment orientations e = (cos(6')) = (AiVs ' Ai)- large L = Ng we have 


1 

P 


- lim - ln((Aiv, ■ Ai)) 

L^OO L/ 


— lim 


1 


Ns —>-oo A Nc 


ln((T 


Ns-l\ 


--i„W 


For K 1, it holds that a ~ /t/3 + O {h?) and hence 




X 

ln(3//t) ’ 


(39) 


(40) 


whereas for k 3> 1 we hnd 

P ^ Xk. (41) 

Therefore, the persistence length is approximated as segment length A multiplied by con¬ 
centration parameter k. 


M2.5 Motor velocity field 

From each realization of the MT cytoskeleton in the polymer model, we derive a nondi- 
mensional vector held termed motor-velocity held v!^ by computing the local vectorial sum 
of MT segments. To this end, we dehne a coarse-grained cubic grid with nondimensional 
side length dG = 0.04 ranging from —1 -|- dG/2 to 1 — dG/2 in x- and y-direction, and 
from Zq -|- dG/2 to Zq — dG/2 in z-direction. The center point for each MT segment is 
computed, and the orientations for all MT segments whose center points fall within the 
same grid volume are vectorially summed to give the average local MT orientation. By 
normalizing the resulting vector held to a mean vector magnitude of 1, we hnally construct 
the motor-velocity held v'^. is used subsequently to compute the cytoplasmic how held, 
and to simulate active motor-driven cargo transport on the cytoskeleton. 


M2.6 Parameter values 

For all computations of the MT cytoskeleton, we used Na = 25000 anterior seeding 
points, giving a total number of points between N = 55764 for wild-type cytoskeleton and 
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N = 84953 for par-1 null mutants. 


The length of an individual MT segment A needs to be short compared to the system 
length L, yet long enough to allow simulations of even long MTs with feasible computational 
effort. Here, we chose A = 0.015 /im. The maximum number of segments per MT is set to 
j\j-max _ 2 gg_ Care was taken that few MTs reach the maximum length A 

Experimental measurements of MT persistence lengths P are on the order of millimeters 
|54j . However, the effective persistence lengths of MTs in oocytes is likely shorter due to 
effects of cytoplasmic flows as well as high density of yolk granules and other obstacles 
in the cytoplasm that induce MT bending. MTs in the oocyte are neither seen to be 
completely straight (Fig. 8, right), nor strongly curled up (Fig. 8, left). Therefore, we chose 


an intermediate value of k = 18, corresponding to an effective persistence length (eq. (41)) 
of P ~ 13.5 /im (Fig. 8, center). 

Due to their curved nature, the reach of a MT polymers is effectively shorter than the 
length of a straight rod composed of the same number of segments. For example, for the 
chosen values of k and A, the end-to-end distance of a MT polymer with N = 25 segments 
(eq. d^) on average 54% as long as a fully extended polymer or straight rod. Therefore, to 
compare the MT polymer model to a model in which MTs are represented as straight rods 
necessitates an adjustment in the mean MT length. 


M3 Rod model for the MT cytoskeleton 

M3.1 Model setup 

In addition to the polymer model for MTs, we test the net orientation of the cytoskeleton 
by comparing with a second model in which MTs are nucleated continuously around the cell 
membrane and act as straight rods with a given length distribution. A single MT from a 
point cr(p, 0) on the boundary that hits an observation point x in the volume contributes 
an orientation vector pointing from cr to a; to the observation point. Heuristically, this 
contribution is expected to be weighted by three different factors: (i) the probability density 
of MT nucleation ps(cr) at boundary point cr(p, 0), (ii) the probability that the nucleated 
MT rod is oriented in a direction such that it hits the observation point x, and (iii) the 
probability that the MT rod is at least long enough to reach from cr to a; across a distance 


26 



To- = ||a; — cr||. Summing up all the weighted orientational contributions from the entire 
boundary eventually gives the net orientation at a given observation point inside the oocyte. 

We denote the joint probability that a MT rod from the surface element dS around a 
point (T reaches a volume element dV around x by p{x,cr)d'EdV. p{x,cr) is normalized 
over the entire oocyte surface area S and volume V as 



didp(a;, cr) 


'V 


1 . 


The marginal distribution p-£{cr) given by 


(42) 


P-eM = [ AVp(x,u) (43) 

Jv 

dehnes the probability density of MT nucleation on the oocyte surface. Conversely, the 
marginal distribution pv{x) 

pv{x) = dT,p{x,(T) (44) 

ds 

is proportional to the total density of MTs that reach the observation point x from the 
entire boundary. Each individual MT rod contributes an orientation unit vector e^a = 
{x — (j)/\\x — cr\\. Thus, the net MT orientation at observation point x corrected for the 
density of contributing MTs can be calculated as 


o{x) 


Pvix) 


dS e^„p{x,(T) 


(45) 


It is convenient to split the joint probability p{x, cr) into the conditional probability and 
the probability for MT nucleation at the surface 


p{x,cr) = p{x\(t) p^{(t) 


(46) 


For the nucleation of MTs on the oocyte surface pj^{cr), we use the same probability density 
as in the polymer model (eq. d^). To specify the remaining conditional probability p{x\(t), 


hrst note that inserting eq. (46) into eq. (43) imposes the normalization condition 


1 = / did p{x\(t) 

Jv 


( 47 ) 
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Consider a local spherical polar coordinate system (ro-, 9„, 0o-) centered at a given point cr 
on the oocyte surface with the Z(j-axis identical to the locally inwards pointing normal h^j. 
MT rods are assumed to be oriented in any direction into the positive half space Zfj > 0 
with equal probability. To determine the general structure of the conditional probability 
density function p{x\(t), we hrst assume MT rods of cross sectional area A and hxed length 
I corresponding to a length distribution 

4>5{ra) = S{1 - r^) , (48) 


with cumulative distribution function 


4>5(r^) 





— r) 


Ta) ■ 


0(r. - 1 ) 


(49) 


The volume Vc = Al of such a cylindrical MT rod can be approximated as integral over the 
positive half space by 


Vc = A dr^ = A dr„Q{l-r„) 


r-27r 


= A 


d(po 


t/2 


d6*o- sin(6*o 


2 0(/ -r^) 


2nri 


^ 1 = / dC 
Jv 


m-rg) 

2nr‘il 


(50) 


where the last step resulted from division by Al. Comparing this result with eq. (47) and 


using eq. (49), we identify the conditional probability density function for randomly oriented 


MTs rods of hxed length as 


p{x\(T)=Mi^^[l-^&{r„)] . (51) 

Herein, the factor Ms = f/l corresponds to the mean MT length, the term (27rr^)“^ rep¬ 
resents the uniform angular distribution on a half sphere, and the expression [1 — <h(ro-)] is 
the probability that MT rods are at least long. 


To account for the experimentally measured MT length distribution, eq. 


(51) is dually 


generalized by replacing ^^(ro-) by the cumulative distribution used in the polymer model 









^T{j'a\n = 3, A) (eq. (31)) with A/r = 1/(2A). Thus, the complete joint probability dis¬ 


tribution p{x,(t) (eq. (46)) for a MT emanating from cr(p) and contributing to the net 


orientation at observation point x in the oocyte volume is dehned as 


p(x,cr) = ps(s\h„, k) Mr [1 - ®r(i'„|3.A)] 


(52) 


Note that MT straight rods always cover the distance between nucleation and observation 
point in a straight line. If the anterior surface is curved inwards, these straight lines can 
temporarily pass outside of the oocyte geometry before reaching the observation point inside 
the volume again. The rod model does not prevent such unphysical contributions (though 
the polymer model does prevent this). However, this error does not occur in a geometry with 
a flat anterior boundary. Calculation of the cytoskeleton in a geometry with flat anterior 
(Fig. 1 - hgure supplement IG) does not change the observed topology, suggesting that in 
practice this error is negligible. 


M3.2 Parameter spaee and eomparison to polymer model 


The continuum description of the rod model allows one to follow the saddle-node bifur¬ 
cation in the MT cytoskeleton precisely throughout parameter space. The creation of stable 
and unstable hxed points always occurs along the anterior-posterior (AP) axis because the 
system is rotationally symmetric. We therefore compute the uncorrected net orientation 


o{x)pv{x) (eq. (45)) for observation points x only along the symmetry axis as a function 


of cortical nucleation parameter kp and mean MT length e. 

For a cytoskeletal topology without any hxed point, the a:-component Ox{x)pv{x) will 
be pointing towards the posterior pole everywhere along the AP axis. A positive value 
Ox{x)pv{x) > 0 everwhere along the central axis therefore dehnes the regions of a wild-type 
cytoskeletal topology. If a pair of a stable and a saddle-node is present, however, Ox{x)pv{x) 
will be negative at at least one point along the AP axis, while at the same time it will be 
positive at the grid point closest to the posterior pole xp where the unstable node creates 
a limited domain of attraction. Combination of both criteria therefore dehnes the region of 
parameter space in which the oocyte is split into two diherent regions of attraction. Finally, 
a negative value at the posterior pole Ox{xp)pv{xp) < 0 is the signature of a cytoskeleton 
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with only a stable point whose domain of attraction covers the entire volume (see Fig. 5). 

In the ensemble averages of the polymer model, we hnd that the bifurcation point occurs 
for small values of /cp if e is large, but shifts to larger values of /cp if e becomes smaller. 
Furthermore, the bifurcation point can be crossed solely by varying e while keeping kp hxed 
(Fig. 4 - hgure supplement 1). In summary, this conhrms that the polymer model exhibits 
the same parameter space structure as the rod model (Fig. 5), thus further supporting the 
similarity of the two models. Moreover, as in the rod model, this parameter space structure 
is again insensitive to the oocyte shape. 

M3.3 Bifurcation normal form 

The saddle-node bifurcation normal form in two dimensions is given by 



wherein A is the bifurcation parameter. For negative A, stable and unstable hxed points 
exist and are located at x± = ±V—A. The critical bifurcation point occurs at A = 0 (Fig. 4 
- hgure supplement IP-T). 

M4 Additional Methods 

Mf.l Nondimensionalization and mathematical methods 

In dimensional units, the system of transport equations for bound and unbound cargo 
fractions c;, and c„ is dehned as 

dtCb + V ■ {vm Cb) = hcu- ku Cb (54) 

4“ • (u Cifj kb Cu T ku Cb ~\~ D ^7 Cu , 

while the Navier-Stokes-Equations and volume forces are dehned by 

pidtu + {u ■ V) u) = —Vp -I- p -I- / 

f Cl Vux , 
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(55) 

(56) 



with no-slip boundary conditions on the oocyte surface. Note that cytoplasmic streaming is 
present in oocytes even when oskar mRNA is not expressed | 32 ], showing that some cargo 
other than oskar mRNA is transported by kinesin and responsible for driving flow. This 


justifies that the parameter a in eq. (56) is independent of the concentration of the bound 
cargo Cb- 

Computations of the fluid flow field and simulations of cargo transport are performed 


using a nondimensionalized version of eqs. (54 56). For nondimensionalization, we select a 


length scale L = 50pm for scaling of space x = Lx', and an advection scale derived from 
the active motor-driven transport V = 0.5 pm/s [H] for scaling of velocities Vm = V v'^, 
implying an advection time scale r = LjV = 100 s for scaling of time t = rt'. v'^ denotes the 
coarse-grained motor-velocity field computed from each realization of the MT polymer model 
(Sec. S2.4). Defining the mean raction rate constant K = {kb + k^) /2 and a nondimensional 


reaction parameter (3 = kb/{2K) we find the nondimensional version of eqs. (54) 


5,.4 + V' ■ (t.;„ 4) = 2Da (Pc', - (1 - ,3)4) 
4-4 + V' . (u' 4) = 2Da (-,34 + (1 - ,3)4) 

+Pe-"V'V„, 


(57) 


where Da = LK/V is the Damkohler number and Pe = LR/D is the Peclet number. 


For nondimensionalization of eqs. (55 56), we define scales for the pressure P = joV/L, 
force density F = pR/L^ and force-velocity scaling A = F/R to obtain 


Re (a; u' + {u' ■ V') u') = - vy + v'y' + f 

/ / / / 

= “ 


(58) 

(59) 


wherein p' = p/F, f = //F, a' = a/A, and Re = pRL/p is the Reynolds number. Note 
that the Reynolds number is defined using the scale R of the active transport field because 
the flow velocity field u is derived from it and hence varies. To estimate an upper bound 
on the Reynolds number, we use the viscosity of water p = 10“^ Pa s and ten times the 
density of water p = lO^kg/m^ as lower and upper bounds for the viscosity and density of 
the cytoplasm in the oocyte. This results in an upper bound Re = 2.5 x 10“^, thus justifying 


to neglect the inertia terms in eq. (58) 
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M4.2 Parameter values 


The unbinding constant ku can be estimated as ratio of the mean active transport velocity 
and the mean track length of oskar mRNA in stage 9 oocytes as ku = 0.17 [H]. The 

nondimensional parameter (3 determines the fraction of cargo that resides in the bound state. 
Given the fact that 13% of oskar mRNA are bound at any given time we set /3 = 0.13 [IT] . 
The unbinding constant ku and (3 can be used to calculate k^ = 0.0255 Therefore, the 
mean reaction rate constant and the Damkoehler number are K = 0.0978 and Da = 9.775, 
respectively. Diffusion of mRNP particles in the ooplasm is widely considered to be very 
slow. Experimentally measured diffusion constants for mRNA transcripts of sizes similar to 
oskar mRNA (2.9 kB) in mammalian cells are in the range of D = 0.02 fim/s or lower [55] . 
With this value, the Peclet number evaluates to Pe = 1250. 

The nondimensional scaling parameter a' (eq. ([5^) which converts the active transport 
velocities into forces acting on the cytoplasmic fluid absorbs several unknown quantities 
including the size of the cargo that drives fluid flow, drag forces on the fluid with local 
viscosity, the density of such cargo-motor complexes on MTs and any collective effects. The 
value of a' determines the mean nondimensional and dimensional fluid flow speeds (|w^|) and 
V (I'u'l). We regard a' as a macroscopic, phenomenological parameter. The value a' = 45 
(a' = 55 for alternative geometry-2) is calibrated such that the resulting mean dimensional 
fluid flow speeds V (|'n'|) match the average fluid flow speeds measured experimentally for 
stage 9 Drosophila oocytes (see Fig. 2) (compare to similar approach in [56]). 

To test the anchoring mechanism, the bound state is replaced by an anchored state to 
which cargo can only bind at the extreme posterior pole. For the anchoring state we set 
ku^^^ = 0 s“^ and = 10 kb, resulting in = 1 (no unbinding), = 0.1275, and 

Da^nch ^ ^2.75. 


M 4.3 Computational methods 

For computation of the cytoplasmic flow field and subsequent transport simulations, the 
motor-velocity held v'^ and forces /' from the MT polymer model are shifted and rotated 
such that the AP axis aligns with the a:-axis and the origin of the coordinate system is 
located in the oocyte center. To compute the cytoplasmic how held we use the open source 
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FreeFem++ finite element solver m- Forces / defined on a 3D cubic grid are interpolated 
to the unstructured 3D grid used by FreeFem++, and resulting fiow velocity vectors are 
interpolated back to the cubic grid and derived staggered grids. 


Cargo simulations of eqs. (57) are performed using custom written Matlab code in finite 


volume formulation on staggered grids as described in |5H] ■ As nondimensional time step we 
used At = 0.005. Each simulation cycles through many pairs of active transport fields v'^ 
and their corresponding cytoplasmic fiow fields u' to account for a dynamically remodelling 
cytoskeleton and temporally changing cytoplasmic fiow fields. Each pair of u^-it'-fields is 
active for 432 time steps, corresponding to 3.6 minutes of simulated real time. 


M4-4 Fly stocks and experimental methods 

Oocyte microtubules were stained as described by Theurkauf et al. [TS] using a FITC- 
coupled anti-alpha-tubulin antibody at 1:200 (Sigma). The general shape and size of the 
oocyte was imaged by recording the fluorescence from a par-1 protein trap (FT) GFP line 
[5^ . Homozygous or heterozygous par-1 FT flies were also used for measurements of cyto¬ 
plasmic flows by imaging movements of autofiuorescent yolk granules. The allelic combina¬ 
tion par-1^^ /par-1^^'^^ was used for imaging cytoplasmic flows in strong par-1 hypomorphs 
[TT] . In all cases of fiow measurements, flies were yeast fed overnight, dissected under Voltalef 
lOS oil and imaged at 40x magnification in a confocal microscope. Cytoplasmic flows were 
recorded with a 405 nm laser, 4 /is pixel dwell time, collecting 13 frames in 25s intervals for 
a total movie duration of 5.17 minutes. 

For the test of oskar mRNA localization under wild-type conditions and under sub- 
physiological temperatures we used an oskMS, MS2-GFP line [HO]- For cold experiments, 
flies were kept at room temperature (22C), yeast fed overnight at 25C, and subsequently kept 
at 25C, 7C and 4C degrees for 6 hours before dissection and fixation. In situ hybridizations 
and imaging of living tissue were performed according to standard methods [32] . The bicoid 
mRNA probe was labeled with Digoxigenin-UTP (Boehringer Mannheim). 
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M4-5 Data analysis 


To compare topologies of in-vivo cytoplasmic flows with topologies of numerically com¬ 
puted cytoplasmic flow fields, we used particle image velocimetry (PIV) to generate vector 
fields capturing the instantaneous direction and magnitude of particle movements. PIV was 
performed using the open source code PIVlab [61], and all PIV vector fields obtained from 
pairs of consecutive movie frames were averaged over the movie duration. For each movie 
of cytoplasmic streaming, the mean flow speeds were calculated as the average over all PIV 
vector magnitudes. Flow speeds obtained from PIV were conhrmed by automatic particle 
tracking using open source code by Blair and Dufresne [ 62 ] • Only particles that could be 
tracked in at least three consecutive frames were accepted. Note that instantaneous speeds 
were computed as ratio of particle distance traveled between frames and frame interval. 
This method avoids underestimation of flow velocities as was pointed out by authors of [26] , 
thereby ensuring that our conclusions about the neglibible contribution of flows to transport 
remain sound. 


M5 Flow fields and oocyte nucleus 

M5.1 Autocorrelation function 

Cytoplasmic flow patterns not only show a high variability between different oocytes, 
but also change over time in any individual oocyte. To estimate the rate of change of 
individual flow helds we calculate the unbiased, discrete vector autocorrelation over time for 
the sequence of PIV flow helds in each movie according to the expression 



(60) 


where angular brackets denote the average across all points of the PIV grid, and the peak 
value of C{k) is normalized to one. Autocorrelation functions from N = 48 movies exhibit 
a large spread across generally low correlation values, and an exponential ht gives a decay 
time constant of 4.4 minutes. 
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M5.2 Nucleus to oocyte volume ratio 


In 2D confocal images, the oocyte nucleus occasionally appears to cover a significant 
fraction of oocytes, particularly in young oocytes up to stage 8. To quantify the size of 
the nucleus compared to the size of the oocyte, we therefore estimate the fraction of the 
oocyte volume that is covered by the nucleus. We first extract the entire boundary of the 
oocyte using a custom written macro in Fiji [63]. Averaging the boundary shape from below 
and above the AP axis results in a symmetrized parametric curve denoted by {x(t),y(t)) 
(Fig. 2 - figure supplement 1C, blue curve). After subdividing the boundary into ante¬ 
rior {xA(t),yA(t)) and posterior {xp(t),yp(t)) curve, the volume of the respective solids of 
revolutions can be computed as 

Vi = Tr J dt {yi{t)f ^ . (61) 


While the anterior surface is typically curved inwards for stage 9 oocytes, the anterior surface 
in young oocytes up to stage 8 can be outwards or inwards curved. Hence, depending on 
the shape of the anterior the volumes of the anterior Va and Vp have to be suitably added 
or subtracted. 

The parabolic caps used in the model geometry (eq. ([^) can be parameterized along the 
2 ;-direction, and in dimensional units their volume is then determined by 


= t: / dz 

Jo 

= ttL^Zq/2 . 



( 62 ) 


Given that the anterior surface in the standard geometry-1 is curved inwards, the total 
oocyte volume is given by the difference between between the posterior and the anterior cap 
volumes Vp —V^ = ttLJ (^Zq — Zq') /2. 

To approximate the volume of the nucleus, we select its shape and fit a circle to the 
selection in Fiji. Using the circle radius, we calculate the volume of the nucleus as 4/3 vrr^. 
We find that the nucleus generally covers less than 2.5% of the oocyte volume (Fig. 2 - 
figure supplement 1), thereby showing that the nucleus volume is negligible even in young 
stage 9 oocytes, and that 2D confocal images convey a poor impression of 3D volumes. 
Given the small volume of the nucleus and its location at the anterior of wild type oocytes. 
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the nucleus is unlikely to influence oskar niRNP transport towards the posterior, and we 
therefore neglect it in simulations of cargo transport. 

M5.3 Impact on flow field 

The nucleus may be expected to disturb the cytoplasmic flow field by means of adding an 
additional no-slip boundary condition inside the volume. We tested this effect by including 
a spherical excluded volume at the oocyte anterior surface in the solution of the 3D Stokes 
equations (Fig. 2 - figure supplement 1A,F,H). The nucleus leads to a clear local disturbance 
of the flow field in a cross section that includes the nucleus compared to the flow field without 
(Fig. 2 - figure supplement 1E,F). However, taking a perpendicular cross section through the 
flow fields in which the nucleus is not visible shows that the nucleus rarely changes the flow 
field (Fig. 2 - figure supplement 1G,H). This low sensitivity with respect to excluded internal 
volumes is due to the fact that flows are driven by volume rather than by boundary forces. 
Note that the flow fields shown in Fig. 2 - figure supplement 1 still likely overestimate impact 
of the nucleus because it generally is not spherical and often seen to squeeze deeply into the 
anterior corners, therefore reaching far less into the oocyte volume than the idealized sphere 
used here. 
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Figure. 1. Models for the MT meshwork show local order in the cytoskeleton. A) 3D geometry 
of a stage 9 Drosophila oocyte (grey, anterior to the left, posterior to the right) containing more 
than 55000 MT seeding points. From the corners to the centre, MT seeding density decreases 
weakly along the anterior {kA = 1000//m, Hq = 0.8) and strongly along the posterior-lateral cortex 
{kp = 150//m, Hq = 0, see Materials and methods M2.1). Nucleated MT polymers are stiff random 
walks, initially pointing in a random direction. Only MT segments in a cross section are shown 
(green) to emulate confocal images. MT target lengths are chosen from a probability distribution 
that accounts for the MT aging process. The mean target length is set to a fraction e of the AP-axis 
length, here e=0.5 (Material and methods M2.3). The inset shows the 3D angular distribution of 
0.5% of all MT segments with 3D statistical bias. B) Cross section through the MT cytoskeleton 
shown in A with 2D directional bias (top right). The inset shows 2D posterior bias (in percent) 
as function of depth (bottom right). C) Local vector sum of MT segments from the cross section 
in panel B on a coarse-grained grid shown as streamlines that visualize local directionality. D) 
Staining of a-tubulin (green) shows MT density distribution in a fixed stage 9 oocyte. Nuclei in 
blue (DAPI), scale bar is 30//m. E) Schematic detailing the work flow in the model and comparisons 
to experiments. F) Local directionality of MT cross section as in panel C for an average over 50 
independent realizations of the cytoskeleton. G) Local directionality computed from the rod model 
with the same parameters as in panels A-B but shortened MT lengths e = 0.25 (Material and 
methods M2.5). Orange arrows show the separatrices between subcompartments. H) MT density 
distribution computed from 50 realizations of the polymer model. 
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Figure. 1 - figure supplement 1. Compartmentalization of the MX cytoskeleton is robust to changes 
in oocyte geometry. A) Alternative geometry for a stage 9 oocyte comprised of a posterior parabolic 
cap and an anterior disc. MX segments intersecting a cross section in one realization of the polymer 
model with nondimensional seeding density parameters kp = 150/rm, = 0,kA = 1000/rm, 

Hq = 0.8 and MX length parameter e = 0.5 from more than 55000 seeding points are shown in 
green. Inset shows directionality of 0.5% of all MX segments with posterior bias. B) Cross section 
through 3D MX cytoskeleton shown in panel A with 2D directional bias (top right). Inset shows 
2D posterior bias (in percent) as function of depth (bottom right). C) Local vectorial sum of MX 
segments visualized as a streamplot. One individual realization of the MX cytoskeleton shows poor 
spatial order. D) a-tubulin staining of an early stage 9 oocyte. Scale bar is 25/rm. E) Schematic 
of the work flow in our model. F) Streamplot of the local vectorial sum of 50 realizations of the 
polymer MX cytoskeleton with parameters as in A,B. G) Local net orientation of straight rod MXs 
in the rod model with parameters as in A,B but reduced mean MX target length e = 0.28, showing 
the same compartmentalization of the oocyte as the average in the polymer model. H) MX density 
distribution in the ensemble of 50 realization of the polymer MX cytoskeleton. 
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Figure. 2. Computed cytoplasmic flow fields capture key elements of in vivo flows. A) Streamlines 
(light blue lines) visualize the 3D cytoplasmic flow held computed from the realization of the 
cytoskeleton shown in Fig. lA. The horizontal plane shows a 2D cross-section through the 3D 
held. Anterior to the left, posterior to the right. B) Cross-section through the 3D held shown in 
panel A with arrows indicating how directions and colouring indicating how speeds. C) Confocal 
image of a live stage 9 oocyte. Arrows show the how held computed from PIV of streaming yolk 
granules and averaged over ~ 5 minutes. Scale bar is 25//m. D) Same as A, but showing the 
mean how organization for an average of 100 individual 3D how helds, analogous to the mean 
organization of the MT cytoskeleton in Fig. IF. E) Same as B for the average in panel D. F) Mean 
huid how speeds were obtained by PIV (red, 13.7 ±0.8 nm/s, mean ± sem) and automatic particle 
tracking (orange, 15.3 ± 0.7 nm/s, mean ± sem) from 48 oocytes. Experimentally measured how 
speeds were used to calibrate the forces f in the Stokes equations Eq. (2) such that the computed 
mean speeds in 3D (blue, mean: 14.5 nm/s) or in 2D cross sections (green, mean: 14.8 nm/s) 
match the measured values. The larger spread in experimental how speeds may rehect greater 
variability of motor activity, cytoplasmic composition, geometry or age in vivo. 
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Figure. 2 - figure supplement 1. The oocyte nucleus covers a negligible fraction of the oocyte 
volume and disturbs the flow field only locally. A) Streamlines visualizing the flow field computed 
in the 3D oocyte geometry with additional no-slip boundary condition on a sphere representing 
the nncleus. Same forces as in Fig. 2A-B, main text. The nucleus occupies 2.2% of the oocyte 
volume. Horizontal plane shows a cut through the flow field. B) Young stage 9 oocyte stained 
with DAPI (blue), phalloidin (red) and showing oskar MS2 GFP (green) in the process of reaching 
the posterior pole. Scale bar is 25/im. C) Extracted and rotated shape of the oocyte in panel B 
(red), and symmetrized geometry resulting from averaging the shape above and below the AP-axis 
(blue). Black circle indicates the nucleus. D) Ratio of nucleus volume to oocyte volume computed 
for Ng = 9 early (red) and = 7 mid stage 9 oocytes (blue). Arrowhead marks the oocyte 
shown in panel B. E-F) Cross sections through the 3D flow fields for the same forces as in panel 
A with (F) and without (E) nucleus. G-H) Same as E and F, but vertical cross section that does 
not contain the nucleus. 
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Figure. 3. The model recapitulates oskar and bicoid mRNA transport, implying dominance of 
cytoskeletal transport. A-C) Shown are fixed oocytes with oskar MS2 GFP (green) and stained 
with DAPI (blue) and Phalloidin (red), oskar mRNA forms a central cloud at late stage 8 (A), a 
collimated channel while moving to the posterior (B), and a posterior crescent at stage 9 (C). Scale 
bars are 25/im. D) 3D oocyte showing the initial oskar mRNA cargo distribution. E-F) Cross 
sections through a simulation of oskar mRNA transport with diffusion, motor-transport and flows 
showing the distribution of total cargo Ch + Cu at the indicated time points. No posterior anchor is 
present. Simulations of 3 (6) hours cycle through 50 Vm-u-'^shis in random order once (twice, see 
Fig. 4). For simulations of 1.5 hours, 25 were chosen at random. Compare to exper¬ 

imental observations in panels B and C. C) Same simulation as in D-F, but without cytoplasmic 
flows, showing largely identical localization as in F. H) Simulation of bicoid mRNA transport shows 
the mRNA quickly accumulating at the nearest cortex (arrowheads) when injected in the posterior 
(inset), corresponding to the behavior of naive bicoid. I) Same as in H for bicoid mRNA injection 
at the anterior-dorsal region (inset). J) Same as in H for injection at the anterior middle (inset), 
showing bicoid mRNA localization to the anterior corners (arrowheads). Localization to the ante¬ 
rior depends on sufficient proximity of the injection site as observed for naive bicoid mRNA. K) 
Same simulation for oskar mRNA as in D-F, but with a posterior anchor (arrowhead) and without 
active motor-driven transport. Compare to panels F and C. 
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Figure. 4. The cortical MT seeding density determines sites of mRNA localisations and gives 
rise to a bifurcation in the cytoskeleton. A-E) Simulations of oskar (A-C) and bicoid (D) mRNA 
transport with diffusion, cytoskeletal transport and flow for the cytoskeletal architectures shown in 
panels E-H, and their corresponding flow fields. Eor oskar mRNA, simulations reproduce wild-type 
localization (A, same as Fig. 3F), and partial (B) or complete mislocalisation (C). Initial condition 
as in Fig. 3D. Simulation times were occasionally increased to 6h (B) to rule out transient con¬ 
centration patterns. For bicoid mRNA, simulations capture mislocalisation to both anterior and 
posterior as in gurken/torpedo/cornichon mutants and in strong par-1 hypomorphs (D, arrow¬ 
heads). Time points as indicated. The inset in D shows initial condition. E-H) Average local MT 
orientations in ensembles of 50 realizations of the polymer model for varying MT seeding densities 
along arclength s (see panel E) of the posterior-lateral cortex (insets). A MT seeding density that 
increases from a wild-type gradient (E, = 0,kp = 150/um) laterally towards the posterior (F, 

Hq = 0,kp = 17.5/im) to a near-uniform distribution (G, Hq = 0,kp = l/rm) shows a saddle-node 
bifurcation by creating a pair of stable (green) and unstable (red) hxed points. H) A MT seeding 
density that is slightly lower at the posterior pole {h^ = 0.7, kp = 40/im) produces mean MT 
orientations virtually indistinguishable from uniform seeding density (compare to G). I-K) Vector 
field of the mathematical normal form of the 2D saddle-node bifurcation (Materials and methods 
M3.3). Values of the bifurcation parameter A as indicated. Fixed points are located at positions 
X = ±\/A. L) Fluorescence in-situ hybridization to bicoid mRNA in a strong par-1 hypomorph. 
The mRNA (red) localises around the cortex, with most accumulation at the anterior corners and 
at the posterior pole (compare to panel D, arrowheads). Scale bar is 25^m. 
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Figure. 4 - figure supplement 1. Either the MT seeding density or the MT length can act as 
bifurcation parameter. A-J) Each panel shows the mean local MT orientations for an ensemble of 
50 realizations of the polymer MT cytoskeleton under variations of the seeding density parameter 
kp and MT length parameter e as indicated. Approximate location of the stable fixed point is 
shown in green, with region of attraction in blue, and domain of attraction of the posterior pole 
in red. The point on the AP-axis between red and blue arrow region marks the unstable fixed 
point. Percentages indicate the the ensemble-averaged 3D directional bias. Experimentally, the 
directional bias was measured as 57.97% : 42.03% m- Cytoskeleton in panel I was used as 
wild-type cytoskeleton (Fig. IF, main text). Anterior seeding density is unchanged in all panels 
kA = 1000/xm, Hq = 0.8. K-0) Same as panels A-J but for ensembles with 100 realizations of 
the cytoskeleton to ensure reliable visualization of vector field even when fewer MTs reaching the 
oocyte center. P-T) Vector fields computed from a saddle-node bifurcation normal form. The 
critical bifurcation point is A = 0, and for A < 0, fixed points are located at x± = ±\/—A. 
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Figure. 4 - figure supplement 2. Lateral MT growth produces the clone-adjacent-mislocalization 
phenotype (see main text). A-D) Addition of lateral MTs reproduces the clone-adjacent- 
mislocalization phenotype. A) Cross section through one realization of the wild-type MT cy- 
toskeleton in the polymer model {h^ = 0.8, kA = 1000/im, Hq = 0,kp = 150/im). Additional MTs 
have been added dorsal to the posterior pole (black arrow) to mimick posterior follicle cell clones 
MARCM) that overexpress dystroglycan [JO]. B) Density of MTs for 50 realizations of 
the cytoskeleton shows enrichment on one side of the posterior pole (white arrow). C) Streamplot 
showing the local vectorial orientation averaged over 50 realizations of the cytoskeleton. Note the 
upwards tilt away from the central posterior pole. D) Simulations of oskar mRNA transport by 
diffusion, cytoskeletal transport and corresponding cytoplasmic flows show that cargo is repelled 
from the site of additional MT nucleation on one side of the posterior pole, thereby capturing the 
CAM phenotype [40] . 
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Figure. 5. The parameter space of the rod model shows the relation and interconversion between 
all three distinct cytoskeletal architectures. A) The regions corresponding to wild-type (red), weak 
par-1 hypomorph (blue) and strong par-1 hypomorph topologies (green, arrowhead at far left) 
are shown as a function of the mean MT length e and extent of the posterior seeding density 
kp. The bifurcation line between wild-type and weak par-1 hypomorph topologies can be crossed 
by either changing the seeding density laterally (horizontal dashed arrow), or by shortening the 
MTs (vertical dashed arrow). B-C) Parameter spaces as in A for increasing MT nucleation at the 
posterior pole (B: = 0.03, C : Hq = 0.1). Inset in panel C shows a magnification of the triple 

point at which small changes in parameters can convert each cytoskeletal architecture into any 
other. D) Local net orientations of MT rods for parameter values indicated in the inset of panel 
C (yellow circles). In all panels, parameter values for anterior MT nucleation were kept constant 
{h^ = 0.8, kA = lOOO/rm). 
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Figure. 6. Seeding point density in the polymer model. Shown are N randomly drawn seeding 
points on the posterior cap of the standard oocyte geometry, according to eq. (15) distributed 
either (near-)uniformly on the cap (A, G) or in a parabolic gradient from the posterior pole to the 
anterior corners (F, L). Between uniform and parabolic distribution, the seeding density can vary 
either by laterally reducing the density (B, C-E, parameter kp), or by reducing the density at the 
posterior pole (H, I-K, parameter h^). Total number of points is calculated with a fixed number 
of anterior points Na = 4000 (eq. (20)). For actual computations of wild-type cytoskeletons 
{kp = 3,hQ = 0, kA = 20, Hq = 0.8, see main text for dimensional parameters), more than 55000 
seeding points are used. 
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Figure. 7. The cortical MT density increases approximately parabolically from the posterior pole 
towards the anterior corners. A) a-tubulin staining of a stage 9 oocyte (green) with DAPI staining 
of DNA (blue). Arrows indicate arclengths sp from the posterior pole to the anterior corners and 
analogously sa at the anterior. B) Fluorescence intensity profile (green) with moving average (red) 
extracted from a 10 pixel wide line along sp indicated in panel A. Arclength was normalized to one, 
and the minimum intensity of the moving average was subtracted from both fluorescence profiles. 
C) Average intensity profiles (red) of cortical MT density as in panel B from = 15 oocytes. The 
minimum intensity value across all profiles was subtracted from each. Blue line shows a parabolic 
fit. D) Probability densities for the distribution of MT seeding points along anterior arclength 
(top) and posterior arclength (bottom) used for the wild-type cytoskeleton (blue lines). Black 
dashed line in bottom panel shows exact parabola. Arclengths are nondimensional lengths with 
scale L = 50/rm. 
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Figure. 8. von Mises-Fisher parameter k determines MT stiffness. Shown are three cross sections 
through MT cytoskeletons for indicated values of k. MTs become stiffer with increasing values of 
K. Insets show angular distribution for 3 x 10^ points drawn from von Mises-Fisher distribution 
around mean direction p, = Bz- 
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